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Abstract 

We perform Hartree-Fock+BCS calculations for even-even nuclei 
with 2 < Z < 82 and ranging from outside the proton drip line to 
the experimental frontier on the neutron-rich side. The ground state 
solutions are obtained for 737 nuclei, together with shape-coexistence 
solutions for 480 nuclei. Our method features the Cartesian-mesh 
representation of single-particle wavefunctions, which is advantageous 
in treating nucleon skins and exotic shapes. The results are compared 
with those of the finite-range droplet model of M0ller et al. as well as 
the experimental values. 

1 INTRODUCTION 



Several theoretical frameworks have been presented for the sake of global 
(on the N-Z plane) calculations of nuclear ground state properties. One 
of the most elaborate results is given by M0ller et al. [|l| using the finite- 
range droplet model with a microscopic shell correction (FRDM). Another 
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extensive calculation is carried out by Aboussir et al. in the extended 
Thomas- Fermi approach. The former as well as the latter methods can be 
regarded as approximations to the Hartree-Fock (HF) equation. While mean- 
field calculations under spherical symmetry are very easily done nowadays, 
those including deformation still requires long computation time for global 
calculations. 

In this paper we show the results of our global HF+BCS calculations 
and compare them with those of FRDM and experiment. Most of the meth- 
ods to treat deformed nuclei [|l], || express single-particle wavefunctions 
by expansion in a harmonic oscillator basis. In contrast, we use a HF-I-BCS 
code ei'(§0, which features a Cartesian-mesh representation of single-particle 
wavefunctions. The mesh representation is expected to be suitable for de- 
scribing the atomic nucleus because the saturation property of nuclear matter 
guarantees the suppression of large-momentum components in wavefunctions. 

The advantage of the mesh representation for the study of unstable nuclei 
is the capability to treat nucleon skins and halos, while they cannot be de- 
scribed correctly in the oscillator basis expansion since the asymptotic form 
of wavefunctions far from the nuclear surface is determined by the basis. 

Another advantage of the mesh representation is that it can treat exotic 
(e.g. high-multipole) shapes and large (e.g. super and hyper) deformations 
without preparing a special basis for each shape. 

2 METHOD OF CALCULATION 

In the code ev8, an octant of a nucleus is placed in a corner of a box (13 x 
13 X 14 fm^), imposing a symmetry with respect to refiections in x-y, y-z, 
and z-x planes (the point group D2h)- The mesh size takes on 1 fm. The 
correction for the inaccuracy due to this finite mesh size is necessary only 
for the total binding energy 0]. (We need a precision of 0.1 MeV in the total 
kinetic energy of 4 GeV for ^°^Pb.) The correction to be added is, 

A^(Z,iV) = 0.00595- (Z + iV) +0.1090- (iV-Z) + 1.24 [MeV], (1) 

which was obtained by fitting to the outputs of a spherical HF+BCS code[|l 
with sufficiently small grid size (0.2 fm). In Fig. |1], the nucleon density of 
a nucleus ^^^Dy in the x-z plane is shown as an illustration of our mesh 
representation. 
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The code uses the Skyrme interaction. Among the widely used parame- 
ters of the interaction, we choose the SIII0. Its vahdity has been examined 
in many nuclear structure calculations. In particular, its single-particle spec- 
trum is in good agreement with experiment. It also reproduces fairly well the 
N — Z dependence of the binding energy compared with other widely-used 
parameter sets0]. On the other hand, its incompressibility is said to be too 
large. 

The pairing correlation plays an important role in determining the nu- 
cleon drip lines It also influences the deformation, sometimes strongly. 
For example, one can see in ref. that the potential energy curve of ^^^Pb 
is substantially changed when the pairing interaction strength is slightly re- 
duced. Therefore, it should be kept in mind that more accurate treatment 
of the pairing is necessary for deformed and/or near-drip-line nuclei than for 
spherical stable nuclei. 

The BCS method cannot be applied to nuclei near the neutron drip 
line because the Fermi energy is so high that positive energy HF orbitals, 
whose wavefunctions are extended over the box, are occupied through pair- 
scattering processes. This brings about an unphysical situation in which 
neutron gas surrounds the nucleus. This problem is resolved in the Hartree- 
Fock-Bogolyubov (HFB) theory, in which ordinary and pair densities are 
localized when the Fermi level is negatively. 

On the other hand, for nuclei near the proton drip line, the HF+BCS 
method is still useful to obtain localized solutions because the Coulomb bar- 
rier keeps the wavefunctions of positive energy HF orbitals localized spa- 
tially. In practice, to prevent the tunneling through the barrier, we modify 
the Coulomb potential outside the barrier so that it is higher than the cut-off 
energy (Ajjp + 7.3MeV) of the pairing interaction. 

We use a seniority pairing force, whose pair-scattering matrix elements 
are defined as a constant multiplied by cut-off factors depending on the single- 
particle energy ef. 

(nKlirln) = -GrfM)fr{e,), (2) 
m = |l + expi^^^^}"VHP + r.3MeV-e),(3) 

where r means proton or neutron while AJjp is the Fermi level of the HF 
(normal) state. For neutrons, the righthand side of eq. is multiplied 
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furthermore by 0{—e). We need a prescription to determine the strength Gr 
for each nucleus in the entire region of the nuclear chart. For this purpose, 
we develop a method based on the continuous spectrum approximation using 
the Thomas-Fermi single-particle level density. As the average pairing gap, 
we adopt the classical empirical formula A = 12 MeV/^/A. 



The imaginary-time time-dependent-Hartree-Fock evolution method[10 



is used to obtain the solutions to the HF+BCS equation. Usually, 500-2000 
steps are necessary to obtain a solution for each nucleus, where the time step 
is chosen as 1.5 x 10~^^ sec. 

In order to search a prolate (an oblate) solution, we first exert an external 
potential proportional to Qz on the initial wavefunction until its quadrupole 
deformation parameter, 

52 = 3(a)/4(r2), (4) 

satisfies ^2 > 0.1 {62 < —0.1). Then, we switch off the external potential, let 
the wavefunction evolve by itself, and see if it converges to a deformed local 
minimum. 



3 RESULTS 

We have performed the HF-I-BCS calculations with the Skyrme SIII force 
for 752 even-even nuclei with 2 < Z < 82 and, for each Z, N ranging from 
several units less than the proton drip line to the heaviest isotope observed. 
The extension of the calculation to the neutron drip line will be deferred 
until a deformed HFB method is developed. Spatially localized solutions are 
obtained for 737 nuclei. Note that some of them are outside the proton drip 
line: The localization is due to the lifting of the Coulomb potential outside 
the barrier. We have also searched for prolate and oblate shape-coexistence 
solutions and found the second minimum for 480 nuclei. 

Although our method allows triaxial shapes, all the solutions we obtained 
are axially symmetric: The 1^2,2 as well as the 14^2 deformations are almost 
vanishing, while the 14,4 deformation is larger but very small. It is possible, 
however, that local minima exist at 7 7^ 0° points in the /3-7 plane which are 
separated from the axial path by a potential barrier. 

Among the quantities which can be extracted from our results are the 
binding energy, even-multipole deformation parameters, radius, skin thick- 
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ness (halo is not found), single-particle levels, and pairing gaps. We also have 
the information on the prolate-oblate difference of these quantities. 

In Fig. 0, the quadrupole deformation parameter \62\ is shown for the 
ground states of nuclei with 42 < Z < 58. One can see the development of 
deformation with the increase of the number of valence particles and holes. 
One can also see that oblate minima are competing with prolate ones for not 
very large deformations {\62\ < 0.2). By comparing Fig. |^ with the result of 
FRDM[|^, we find conspicuous differences for 40 < Z,N < 50, where only 
the latter predicts very large deformations. 

In Fig. ^ two-proton separation energies are plotted for nuclei in the 
neighbourhood of ^°°Sn. The results of FRDM|]l| and experimental values [|lT 



are also shown. The agreement between them looks more impressive consid- 
ering that the SIII force was made in 1975 while FRDM was revised in 1993. 
A noticeable discrepancy is that the jump across the Z = 50 major shell 
closure is smaller in experiment than our calculation. 

In Fig. the thickness of the proton skin is shown. The definition of the 



skin is taken from ref. |]T2|: The condition for a point r to be in the proton 
skin is, 

Pp(r) > 4p„(r) and pp(r) > Pp(0)/100. (5) 

When the nucleus is deformed, the thickness of the skin can depend on the 
direction. In the figure the skin thickness in the z-axis (symmetry axis) is 
plotted. The skin in the x-y (equatorial) plane is usually slightly thinner 
for prolate shapes because of the Coulomb repulsion. One can see that the 
proton skin is very thin (<0.3 fm) for nuclei inside the proton drip line except 
some light {Z < 20) nuclei. For all the nuclei with 50 < Z < 82, which are 
not shown in the figure, the proton skin does not exist. 
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FIGURE CAPTIONS 



Fig. 1. The nucleoli density of ^^^Dy in the x-z plane calculated with the 
HF+BCS method with the SIII force. This nucleus is located just in- 
side the proton drip line and has a prolate shape (^2 = 0.315) which is 
symmetric in the 2;-axis. Dashed lines correspond to the 1 fm mesh. 

Fig. 2. The quadrupole deformation parameter \62\ as a function of the neu- 
tron number N for some isotope chains near Sn calculated with the 
HF+BCS method with the SIII force. On the lefthand (righthand) side, 
Z = 42,44,46,48 {Z = 58,56,54,52) isotopes are designated by trian- 
gles connected with solid lines, squares with short-dash lines, diamond 
shapes with long-dash lines, and circles with dot-dash lines, respectively. 
The solid symbols are used for prolate states, while open symbols are 
used for oblate ones. All the Z = 50 isotopes are spherical and have no 
deformed local minima. 

Fig. 3. Two-proton separation energy in the neighborhood of ^'^"Sn. The cal- 
culation with the HF-I-BCS method with the SIII force is designated by 
solid circles. The solid lines connect the solid circles of the same neutron 
number. The neutron number is printed at an end of each line. The re- 
sults of FRDM|l[] are represented by open circles, while the experimental 
values are designated by plus marks. 

Fig. 4. The thickness of the proton skin in the 2;-axis (the symmetry axis) 
calculated with the HF+BCS method with the SIII force. The definition 
of the skin is given in the text. The solid (open) circles are used for nuclei 
with S2p > (S'2p < 0). The lines connect isotones. The number printed 
near each line means the neutron number of the isotone chain. 
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